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During inflation, scalar fields, including the Higgs boson, may acquire a nonzero vacuum expecta¬ 
tion value, which must later relax to the equilibrium value during reheating. In the presence of the 
time-dependent condensate, the vacuum state can evolve into a state with a nonzero particle num¬ 
ber. We show that, in the presence of lepton number violation in the neutrino sector, the particle 
production can explain the observed matter-antimatter asymmetry of the universe. We find that 
this form of leptogenesis is particularly effective when the Higgs condensate decays rapidly and at 
low reheat scale. As part of the calculation, we present some exact results for the Bogoliubov trans¬ 
formations for Majorana fermions with a nonzero time-dependent chemical potential, in addition to 
a time-dependent mass. 


I. INTRODUCTION 

During the inflationary era, the Higgs field may de¬ 
velop a stochastic distribution of vacuum expectation val¬ 
ues (VEVs) due to the flatness of its potential [IH3]) or 
it may be trapped in a quasi-stable minimum. In both 
cases, after inflation the Higgs field relaxes to its vacuum 
state via a coherent motion [3HS]- At large VEVs, the 
Higgs field may be sensitive to physics beyond the Stan¬ 
dard Model, and new terms in the Lagrangian, such as 
those considered in models of spontaneous baryogenesis, 
can generate an effective potential for baryon and lepton 
number [3IS]. These terms couple the time-dependent 
scalar condensate to the lepton (and baryon) number cur¬ 
rents. Consequently, the thermal bath of particles pro¬ 
duced by reheating can lower its energy by converting 
particles into antiparticles, through scattering involving 
neutrinos, whose Majorana mass violates lepton num¬ 
ber [51 [5]. This scenario can explain the observed matter- 
antimatter asymmetry of the universe; it requires suffi¬ 
ciently fast reheating, such that the plasma forms before 
Higgs relaxation is complete, which restricts the possible 
parameter space. In addition to the Higgs boson, an ax- 
ion or a Majoron relaxation could generate the baryon 
asymmetry of the universe liiaiiQ]. 

However, the relaxation of the Higgs vacuum expecta¬ 
tion value itself results in particle production. Generi- 
cally, a classically evolving background scalar field cou¬ 
pled to quantum fields results in particle production; this 
can be understood as a consequence of the fact that the 
initial vacuum state (which is annihilated by the appro¬ 
priate annihilation operators at t = 0) is not annihi¬ 
lated by the appropriate annihilation operators at later 
times [TTI - H3] . (More specifically, the time-dependent 
background mixes positive and negative energy solutions 
of the field equations, and so an initially diagonal Hamil¬ 
tonian is non-diagonal at later times.) This can result 
in the production of both scalar bosons O [13] and 
fermions miii], provided that the classical scalar field is 


coupled to both. This has been explored extensively with 
respect to the inflaton (e.g., [T5lH8) L In this work, we 
calculate the excess of neutrinos over antineutrinos pro¬ 
duced by the evolving Higgs condensate in the presence 
of chemical potential, generated by higher order terms 
in the Lagrangian, which distinguishes particles from an¬ 
tiparticles. 

During the oscillations of the Higgs condensate, the ef¬ 
fective chemical potential changes sign, which alternates 
whether the production of neutrinos or antineutrinos is 
favored. Therefore, the maximal asymmetry is produced 
with the Higgs condensate decays quickly, which mini¬ 
mizes this wash out. Furthermore, this mechanism favors 
a low reheating scale, which minimizes entropy produc¬ 
tion. 

The outline of this paper is as follows: In the next 
section, we introduce our model, including the Oq opera¬ 
tor which gives rise to an effective chemical potential for 
lepton number. Subsequently, we derive an effective La¬ 
grangian by integrating out the weakly interacting right- 
handed neutrino states, and we specialize to the case of 
a single fermion family. In section |IV[ we quantize this 
system and find the Bogoliubov transformation equations 
which describe particle production. Following this, we 
define the occupation number of the physical eigenstates 
and lepton number. Finally, present a numerical analysis 
of our model, which demonstrates that resulting asym¬ 
metry can be sufficiently large to account for the observed 
matter-antimatter asymmetry. 


II. THE LAGRANGIAN 

In this section, we introduce the model Lagrangian. 
We begin with the action in general curved spacetime 

S' = y d'^x^/^C 


( 1 ) 
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with a Lagrangian 

= £//+/Cf + + £sm, ( 2 ) 

where we use Ch to denote the Higgs sector contribu¬ 
tion, Cf, to denote the lepton sector contribution, /loe 
to denote higher dimensional operators which will gener¬ 
ate an effective chemical potential for baryon and lepton 
number, and £sm represents the Standard Model contri¬ 
butions that do not appear in Ch or We consider an 
expanding FLRW spacetime with signature (-I-, —, —, —). 

The purely Higgs sector contribution is 

CH=9^''d^^^d,^-V^{<^,T), ( 3 ) 

where is the Higgs potential, including any rel¬ 

evant loop and finite temperature corrections. We note 
that, as with the models discussed in mm, the poten¬ 
tial may require higher dimensional operators 

involving the Higgs field d) (and possibly the inflaton 
field I) in order to suppress isocurvature perturbations 
resulting from variations in the produced baryon den¬ 
sity [TM^ . The Higgs sector is discussed in more detail 
in subsection III Al below. 

The lepton sector Lagrangian includes the terms 

Ci = i^L -f L 

+ i^NR “ X! 

- Y, V.L^Nr - Y ^WrTNr + h.c., (4) 

where L denotes left-handed lepton doublets, £ right- 
handed charged leptons, and Nr right-handed neutrinos, 
and we implicitly sum over indices and families. The 
gamma matrices in the FRW metric are related to those 
in flat space time by 7 ^ = 07 ^. We note in Eq. @> 
the effect of the spin connection evaluated on the FLRW 
background. This Lagrangian will be discussed further 
in subsection III HI below. 

The third part of the Lagrangian is the higher dimen¬ 
sional operator 


^Oe — (5) 

where is the lepton current density. One possibility 
for generating an operator of this form is to couple the 
Higgs field $ to the SUl(2) x Uy( 1) gauge fields A and 
Bby 

( 6 ) 

which can be written in the form of using the electroweak 
anomaly equation mm- This transformation requires 
the electroweak sphalerons to be in thermal equilibrium, 
which may not be satisfied here, although the situation 
is complicated by the time-dependent Higgs VEV. For 


these reasons, we discuss other ways of generating this 
operator in Appendix [Aj 

We will discuss the role of this term further in subsec¬ 
tion m for now, we note only that in the presence of 
this term, the Higgs evolution induces a chemical poten¬ 
tial that distinguishes particles from antiparticles. 

Before we discuss each component separately, we will 
first rewrite the action using conformal time, such that 
the metric is = 0 ^( 77 ) 77 ^^, ancQ 


/■* dt 

(7) 


It will be convenient to define the “comoving” 

(j) = acj) 

fields. 

-ijj = 

such that we can write 

( 8 ) 

S' = y d!^x (^Ch Ci + Coq + 

(9) 


where 
Ch = 

Ci = iYL0L + iY Nr$Nr - Y y'^^^R 

- Y vL^Nr - Y ^{NrYNr + h.c., 

•^SM = (10) 

In the first equation, we have defined a comoving poten¬ 
tial Vtf, = a^V^. 

In the next subsection, we consider how the Higgs field 
might acquire a large VEV during inflation, which re¬ 
laxes to its equilibrium value during reheating. Then 
we consider how this affects the quadratic terms in the 
lepton sector; subsequently, we demonstrate that when 
the Higgs VEV is in motion the Og operator produces 
a chemical potential for baryon and lepton number. Fi¬ 
nally, we gather together the relevant contributions to 
the Lagrangian in the final subsection. 


A. The Higgs Sector 

The Standard Model Higgs boson has the tree-level 
potential 

V^($) (II) 


^ Throughout this paper, we will use primes to denote differentia¬ 
tion with respect to conformal time, dots to denote differentiation 
with respect to physical time, tildes to denote comoving quanti¬ 
ties, and hats to denote two-component fields. Where necessary, 
we will use bars in dummy variables. 
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where $ is the Higgs SU(2) doublet. The parameters 
m and A, although constant at tree-level, are modified 
by both loop and finite temperature corrections. For 
the experimentally preferred top quark mass and Higgs 
boson mass, loop corrections result in a negative run¬ 
ning coupling A at sufficiently large vacuum expecta¬ 
tion values (VEVs), with the result that the = 

■yEM = 246 GeV minimum is metastable at zero temper¬ 
ature [22. We note, though, that a stable vacuum is pos¬ 
sible within current experimental uncertainties [22] . and 
the stability of the potential is also sensitive to Planck- 
scale corrections [23] ■ 

Therefore, the running quartic coupling produces a 
shallow potential, and consequently, the Higgs field may 
develop a large VEV during inflation due to quantum 
fluctuations |4] . Qualitatively, the scalar field in a de Sit¬ 
ter space can develop a large VEV via quantum effects, 
such as Hawking-Moss instantons [D [23 or stochastic 
growth [2123123. Subsequently, the field would relax to 
its equilibrium value via a classical motion on the time 
scale ~ unless Hubble friction delays this 

relaxation. If Hj ^ y/d?Vjd(j? then quantum jumps oc¬ 
cur frequently enough to maintain a large VEV. 

Alternatively, the Higgs potential is sensitive to higher¬ 
dimensional operators at large vacuum expectation val¬ 
ues, which can have the effect of lifting the second min¬ 
imum, stabilizing the electroweak vacuum. During infla¬ 
tion, the Higgs field may have a stochastic distribution of 
VEVs similar to that of the inflaton itself in chaotic in¬ 
flation models. During inflation, sufficiently large VEVs 
evolve towards the false vacuum from above, and then 
remain trapped in this false vacuum until destabilized 
by thermal corrections during reheating. Subsequently, 
the field rolls to the global minimum, until electroweak 
symmetry is broken at a significantly later time. 

Therefore, it is quite natural to consider scenarios in 
which the Higgs field has a large vacuum expectation 
value during inflation, which subsequently relaxes to its 
equilibrium value. Both of the above scenarios have been 
explored previously [3I3- Here, we consider the classical 
motion of the Higgs field towards equilibrium generically, 
without specifying the mechanism which generates the 
initial large vacuum expectation value. 

We note that if the field has expectation value 


($) 


1 

71 



( 12 ) 


then the comoving field has expectation value 


$ 


72 V 0 y 


J_ {av{T])\ 

72I 0 J 


(13) 


where we have defined v = av. For completeness, we dis¬ 
cuss the equation of motion for v in the Standard Model, 
including loop and finite temperature corrections, in Ap¬ 
pendix]^ 


B. The Neutrino Sector 


Next, we consider the effect of the evolving Higgs VEV 
on the quadratic terms in Ci, given by the second line of 
equations (10). Including multiple generations, we write 
this as 


Ci = i ^ La$La + i ^ Njii0Nfii — E yia^L 

aa 

ot i a/S 

— yvajeabLga^bNRj — + h.C., 

aj ij 

(14) 


where $ is the comoving Higgs doublet, L is the comov¬ 
ing left-handed {vlYl) lepton SUl doublet of species L, 
and N are right-handed Majorana neutrino states. Greek 
indices label flavors (e, /i, r), while the Latin indices i 
and j label right-handed neutrinos. The indices a and b 
are SUl labels. These are the only renormalizable terms 
which describe the interactions between the Higgs and 
lepton doublets, given the gauge symmetries of the Stan¬ 
dard Model. 

When the comoving Higgs field acquires a vacuum ex¬ 
pectation value, this becomes: 


£g = i'^ La$La + i - 7^ 


yea/3V 


^Lol^R^ 


-E 


yuaj^— - 


72 


-vl, 


t^Ri — 7 ^ 


(NmYNjij + h.C., 


(15) 


where the comoving mass is M = aM. The right-handed 
Majorana neutrinos induce lepton-number violation in 
interaction involving neutrinos; however, there is no cor¬ 
responding effect for the charged leptons. Therefore, 
these terms will not affect our analysis, and so we will 
absorb them into £sm- We define the neutrino sector 
Lagrangian 

Cu=i'^ VLa$VLa + * Nr^Nih - Mj^NRiVLa 

OL i Oil 

— 2 E ^NijiNRiYNiij + h.C., (16) 

ij 

where the comoving Dirac mass is 

(17) 

We note that since v = av, this has the expected scaling 
of a comoving mass. 

It will be convenient to use two-component comoving 
Weyl spinors; we work in the chiral basis, with conven¬ 
tions outlined in Appendix [C| 
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We denote the two component spinors with tildes, as 

VL = (4 0 ) = (0 vl) , 

^=(0 ^i)(? J) = (^^ 0). (18) 


by parts by moving the derivative onto the Higgs vacuum 
expectation value. As we are in flat FLRW spacetime, we 
may replace with to find 

^Oe = ^2 ^B+L- (25) 

Finally, we want to express this in terms of the comoving 
fields. We note that the current is 

3^iB+L = (26) 


The conjugated comoving fields are 

Ji = CN~r = (0 -iNla 2 ) , (19) 

where the charge conjugation operator C is also given in 
Appendix 

In terms of these two-component comoving spinors, the 
neutrino sector Lagrangian can be written as 

^ + i X] 

ol i 

Oii 

( 20 ) 

and we note that if we introduce the fields 

Ni = -ia2N*R ( 21 ) 

the mass term has the expected Majorana form, 


where we have defined jfiB+L = intro¬ 

duce the notation of a prime to denote a derivative with 
respect to conformal time. This allows us to write 


where 



2a'$^Jo^/a -^ 

M 2 ^b+l- 


= 2aa'$^ + 


0/.A 




(27) 


(28) 


When $ acquires a time-dependent vacuum expecta¬ 
tion value u, this is 


^Oe — - 


dov 


where we emphasize that 


^ — 2a'v^ja~Q 
2M2 ^b+l- 

do = d/dt]. We define 


(29) 


= 


dov"^ — 2a!v^ ja 
' 210^ ' 


(30) 


since this term acts like a chemical potential for H -|- L 
charge. We note that in terms of v and M, this has the 
expected form 


- + N\Ng), (22) 

where we remind our readers that these are comoving 
fields. 


C. The Oe Operator 

Next, we turn our attention to the Og operator intro¬ 
duced in equation 0; we use the lepton and baryonic 
currents 

7^=E-'?V9, (23) 

q I 

where the sums are over all leptonic fields, including 
right-handed neutrinos, and baryonic fields respectively. 
In a general curved spacetime, we assume this becomes 

to, = (24) 

which generally holds if the gravitational anomaly is can¬ 
celled by having equal numbers of left- and right-handed 
neutrinos, as discussed in Appendix]^ Next, we integrate 


a =-^ —02 _ 2—-02 

2M2 ydr] drja J 

a dv 

^ ■ 

Therefore, we also define 


(31) 


1 dv 


(32) 


such that p, = afi. We note that since the Higgs VEV is 
initially decreasing, dv/dt is initially negative. The neg¬ 
ative sign in Oq was chosen in order to bias the creation 
of particles over antiparticles. 

As mentioned previously, lepton number is violated 
only in the neutrino sector, and therefore we are inter¬ 
ested only in contribution to the current from the neu¬ 
trinos. Thus the relevant part of the Oo operator is 




+ 0 ) 




Oc 


(33) 
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D. Complete Two-Component Neutrino 
Lagrangian 

Using the results of the previous subsections, the 
complete effective Lagrangian for the comoving two- 
component neutrino fields is 

cx i 

- ^ + A Am 

Oil 

- \ (iMAIAr, - A^RAm 

+ A ^ {yLa^La + i (34) 

Ol 

which will be the basis for our subsequent analysis. This 
describes a set of left-handed and right-handed neutrinos, 
with a Dirac mass and a right-handed Majorana mass, 
and a chemical potential for neutrino number, obtained 
from (IMl). 


convenient to integrate out the heavy right-handed neu¬ 
trino^ and specialize to a single generation, which will 
be a sufficiently rich model to capture the asymmetry 
production of interest here. 

The comoving right-handed neutrinos obey the equa¬ 
tions of motion 


0 = ia^d^Nm - ^ MoiaVLa + i cr2^Hi 

a j 

+ jlNRi. (35) 


In the limit of small p, and when the kinetic term is neg¬ 
ligible, namely at scales below that of the right-handed 
Majorana mass eigenvalues, this equation is solved by 


III. THE EFFECTIVE LAGRANGIAN FOR ONE 
GENERATION OF LEFT-HANDED NEUTRINOS 


Njik = —i 


(36) 


The Lagrangian in Eq. (34) includes several genera¬ 


tions of both left and right-handed neutrinos. It will be which when substituted into the Lagrangian gives 


£eff = -'-Y. 

a a,(3 

+ ^D)a 0 i'Ra^Lp- 

Q(,/3 

Doing so induces a Majorana mass for the left-handed neutrinos, of magnitude, 

Ml = MIA^Md. 

Thus this Lagrangian may be written (using implicit notation for the sums): 


Loc 




(37) 


(38) 


■^eff — ^n^La 


2 L 


{ML)a/3i>la<^2C'L0 - {MI)AlAr 


"b ffi ^D)aj3R\ja^LPi 

(39) 

which has the equations of motion, 

0 = ia^d^i>La + i Y(^L)a/3Cr2t>lp + (40) 


^ We do emphasize that we integrate out Nji, which are not, 

strictly speaking, identical with the heavy mass eigenstate. This is a good approximation below the scale of the right-handed Ma- 
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It is beneficial at this point to specialize to the one- 
generation case, since, as we will show, one generation 
is enough to obtain a nonvanishing asymmetry in the 
presence of the Oq operator. We see explicitly that the 
induced Majorana mass transforms with the form appro¬ 
priate to a comoving mass. 


Mr. = 


M 


Mn 


2 2 

D_ ^ V 
2Afjv 


(41) 


Thus we use the effective Lagrangian, 


CeS = iv\a^duVL - 


iMl 


i'1(72VL 




(42) 


where 


Meff — 






(43) 


wheen = ylv'^/2M‘^ ^ 1. We have rotated the 

field to eliminate the phase in Ml which arises from 
the phase in the Yukawa coupling yi,. (Note that the 
Higgs VEV V can be taken to be real at all times.) 


IV. QUANTIZATION AND BOGOLIUBOV 
TRANSFORMATIONS 


Let us now discuss quantization. First, we consider 
the scenario in which the mass and chemical potential 
are time-independent; we solve the equations of motion 
and determine the creation and annihilation operators 
which diagonalize the Hamiltonian. 

Then we include the time-dependence of the mass and 
chemical potential, which induces a mixing between the 
positive and negative energy solutions of the field equa¬ 
tion. Consequently, even if the Hamiltonian is diagonal 
at time t = 0, at a later time it will be non-diagonal. 
It may be diagonalized with a time-dependent redefini¬ 
tion of the creation and annihilation operators; the co¬ 
efficients of this diagonalization are known as the Bo- 
goliubov coefficients, and in the subsequent section, we 
will relate these coefficients to the occupation number 
of physical eigenstates and to the lepton number. This 
follows the procedure of e.g. [STHSS!- 

We do note that in the multi-generation case, the time- 
dependent rotation that diagonalizes the mass matrix can 
introduce novel effects into particle production, as dis¬ 
cussed in [29]; however, such features will not be neces¬ 
sary to generate a nonzero lepton number. Therefore we 
work in the one generation limit, using Eq. (42), where 
these terms are absent. 


jorana mass eigenvalues. 


A. Constant Mass and Chemical Potential 


The equation of motion for Lagrangian with a single 
comoving Weyl field, Eq. (42), with constant comoving 


mass Ml and comoving chemical potential jleS.i is 

{idr) - i(T ■ d)vL = -ML{i(T2)vl - jleSVL- (44) 
or in momentum space, 

(iSo -f h\k\)vL = -ML{i(T2)i'*L - iieei>L- (45) 

For consistency, we will use k for the comoving momen¬ 
tum, and p for the physical momentum. We consider a 
solution of the form 


12L = 


d^k 

(^ 






h=±l 


, (46) 

where is the two-spinor which is an eigenstate of 

the helicity operator (appropriate to k) with eigenvalue 
h = ±1. This ansatz, when substituted into the equation 
of motion, requires 

{idrr h\k\)u{h, k) = hMLv{h, k) — fLesu{h, k) 

{ido + h\k\)v{h, k)* = —hMLu{h, k)* — flegv{h, k)*. 

(47) 

These equations can be decoupled, 

{ido + h\k\ + p,es)i-ido + h\k\ + fles)F = -MIF (48) 
with F = u,v, satisfies 


d^oF = 


{h\k\ + p-es)^ + Ml 


F. 


This has solutions of the form F = where 


w = y {h\k\ + /Ieff)2 -I- Ml- 
Therefore, we take 

u{h, k) = + 


(49) 


(50) 


72 


72 


;(/r,7 = 


(51) 


where 


/= (52) 

Cj 

and a and j3 are constant coefficients. (In the time- 
dependent case, these will be the Bogoliubov coeffi¬ 
cients.) One can verify that these satisfy the first order 
equations of motion. 
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The state vl obeys the anticommutation relations 
{vL{x),VL{y)‘^] = y), 

{i>L(a;),z>L(y)} = 0, ^^i'l{x),i>l{y)^ = 0. (53) 

These follow from the ansatz 


{af),af}=0, {af\af^}=0. 


(54) 


along with the normalization condition |ap + |/3p = 1. 

Next we proceed to diagonalize the Hamiltonian; the 
appropriate creation and annihilation operators will not 
be the and operators themselves, but linear 

K K 

combinations of these operators. Note that even with 
the Oq operator, the equation of motion ensures £ = 0; 
this follows from the fact that the equation of motion is 
first order. Therefore the Hamiltonian is 


d^x (j>ldoi)L - . 

In terms of the a operators, this Hamiltonian is: 


^=2 


(55) 


-4 


d^k 

(^ 


[2 [|a|2 - 




+2ha*(3*C{k, + 2/ia/3C(fe, /i)*a4a4 

(56) 

where we have introduced the notation pu for the four- 
vector {E, —k). C(fc, h) is a phase factor which arises from 
the product of the two spinors. For the interested reader, 
the important steps in this derivation are discussed in 
Appendix This can be written as a matrix equation: 




d?k 

(^ 


^“(“r’ »£) 


( l«l“-|/3|“ 
\2ha(3C{k,hY 



2ha*Y*C{k,h) 


We introduce the rotated states: 


^(A)tN 


a* hYC{k, h) 
-hY*C,{k,h) a 



(57) 


(58) 


which diagonalizes the Hamiltonian 




{2ttY ^ ^ ku ko ' 


d^k 

(27r)3 




(^)t 

k '^k ’ 


(59) 


where we have normal ordered and changed the integra¬ 
tion variable to —p in the second term. We note that 
since h^ = l, we can also write the eigenvalues as 


u}= J (|fc| -I- hpesY + Ml. 


(60) 


Additionally, we note that as expected, Cj = auj^ where 

w = ^J{\p\ +hp^sY + Ml, (61) 

and p is the physical momentum corresponding to the 
comoving momentum k. 


B. Time-Dependent Mass and Chemical Potential 

Now we consider the case in which both the comoving 
mass and comoving chemical potential evolve in time. We 
again use an expansion of the comoving Weyl spinor of 
the form of Eq. ( [4^ , and the equation of motion again re¬ 
quires u and V to satisfy equations of the form of Eq. (471, 
but with time dependent quantities Ml and peg. 

We will consider solutions of the form 


(h, 4 = E 4We-* /c 


V2 

ha 


V2 

hl3 


(h, ~k) = y ^ 


^ Codr} 


^ Codr} 


V2 


V2 


(62) 


As shown in Appendix E this leads to the differential 
equations 


da 

df] 


dY 

dr] 


2 £2 




dr] 


dr] 


pii fg'i cudr/ 


a 1 
2 £2 




• e 


dr] 

-2z/J' Qdfj 


dr] 


(63) 


(64) 


We take the initial conditions to be a{r] = 0) = 1, 
/3{r] = 0) = 0. This is consistent with the normalization 
condition \al + |/?p = 1, and at t = 0, the A operators 
align with the a operators, so the Hamiltonian (at this 
time) is diagonal when expressed in terms of either set. 
The diagonalization of the Hamiltonian proceeds as in 
the time-independent case, as discussed in Appendix]^ 
The effect of the time-dependent comoving mass Ml 
and chemical potential p^s is to mix positive and nega¬ 
tive frequency modes, as is evident by the fact that (3 will 
generally be nonzero at later times. From the transfor¬ 
mation matrix Eq. (58), the operators that diagonalize 


the Hamiltonian at later times will generally be nontriv¬ 
ial linear combinations of and 


V. PARTICLE NUMBER AND LEPTON 
NUMBER OPERATORS 

Next, we express the expectation values of the occupa¬ 
tion number operator (for the physical eigenstates) and 
the lepton number operator in terms of the Bogoliubov 
























coefficients a and /3. As the operators and 
diagonalize the Hamiltonian, these correspond to physi¬ 
cal particles. The procedure that we follow is this: we 
first express Nh and Leff in terms of these operators 
and normal order (for a discussion on normal ordering 
see [30]). We then express the operator in terms of the 


aip'^ and operators using the transformation equa¬ 
tions (58). Then, we take the expectation value with the 
state f\WC, 0), the vacuum at time t = 0. 


The total number of physical particles of helicity h is 


^ (vAC; 0|Af )|VAC; o) . (65) 


This operator is already normal-ordered, so we proceed 
to write this in terms of the time-independent oper¬ 
ators. 


form annihilate the vacuum. Therefore, 

k 

(JVAC; 0| |VAC; 0^ 

= |/3p (vAC;0|a[,'‘^a[;''^VAC;0\ . (67) 

\ fc£) «£) / 

This matrix element is 
(vAC;0|agaWt|VAC;0) 

= (27r)®5(®^(0)<5h,h (VAC;0|VAC;0) - (vAC; 0|4'‘'^4 '*VaC; 0 
= Vcm, (68) 

where Vcm stands for the comoving volume (and we 
have used the usual formal manipulation (27r)^(5(0) = 

/ ° = Vcm)- This gives the expected result 

/ d^k 

(69) 




+ hal3({k, h)'afyf‘’ + 

( 66 ) 


where in general, 13 may depend on the momentum and 
helicity, as we have noted. 

Next we consider effective lepton number, which is car¬ 
ried by the neutrinos. This charge is given by 

Leff = [ d^xvlvL- (70) 


We assume that we are in the state |VAC, 0), the vac¬ 
uum state at time t = 0. Therefore, all operators of the 


Following the procedure outline above gives us a normal 
ordered expression 


Teff — 


d^k 

(^ 


E 


k k ^ ZUJ ^ 


(71) 


where the important steps are described in Appendix]^ Taking the inner product with the t = 0 vacuum gives 

d^k 


VAC;0| : : |VAC;0) = 


( 27 r) 




-/) (VAC; Q\Af^'^Af'^ |VAC; 0 


Ml 

2uj 


2i/j'id^Wy^C;0|A('‘^U('‘^VAC;0') - (vAC; 0|A('‘^Af ^|VAC; 0 

\ k k o j 2uj \ k]j k 


(72) 


We express these in terms of the operators, which an¬ 
nihilate the state |VAC;0); however, now that we have 
normal ordered we are careful to maintain any Dirac delta 
functions that arise from using the anticommutation re¬ 
lations. The second and third matrix elements are 


(vAC;0|Af)^A('^^VAC;0) = ha* K* Vcm 
(vAC; 0| aI'^^ A^ > |VAC; o) = hap* (Vcm (73) 


Therefore, the lepton number as a function of time is 


(vAC; 0| : Zeff : |VAC; o) = Km 


d^k 

(27r)3 


E [(-/)i/ 3 'i 


- — h (a*Pe^^^o^dfi y 

2lu \ J 


( 74 ) 


where / = {h\k\ + pes)/ui. 
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VI. ROTATED OPERATORS 

At this point, it is convenient to define the rotated 
operators: 

a = ae-^ 

P = pe^IS^df) 


(75) 


which obey the differential equations 

^ = —c{ri)p — iuja 
dr] 

dp _ - 

— = c(r])a + iwld, 


where 


1 1 

2^ 


Ti~f dpes / LI r I I ~ ^ 

Ml—, - {h\k\ +Meff)—^ 

dr] drj 


(76) 


(77) 


and we also have the normalization condition + = 

1 along with the initial condition d = 1 and /3 = 0. 

Note that we can rewrite this so that the helicity h 
multiplies the chemical potential, 


h 1 


MLh^-i\k\ + hil,s)'^^^ 


dr] 


dr] 


(78) 


In terms of these rotated coefficients, the number den¬ 
sities and lepton number are: 


VAC; 0| : : |VAC; 0^ = 14,, 

VAC; 0| : Leff : |VAC; o) = 14, 

H- /ieff 
Co 


(Pk 


sE 


(27r)3i^fc./»i ’ 

" d^k 

{2tt) 


Mr 


H,h\-^^C^l,hPk,h + ^k,hPlh 


(79) 


which has eliminated the fast oscillatory time depen¬ 
dence. The comoving number densities are therefore 


nu = 


riL = 


— Id- |2 

^2Trr\Pk,h\ , 

d^k 

(27r)3 


E 


h\k\+lleS.3 |2 Ml 


\Pk,h\^ 24 ; ^ ip^k,h^k,h + ^k,hf^k,h 

(80) 

At late times, the Higgs VEV v approaches zero, and 
therefore the comoving VEV v = av also does. Conse¬ 
quently, Ml —t 0, /ieff —t 0, and Co —)■ |fc|. Therefore, the 
limit of the comoving lepton asymmetry is 


lim fiL = 


d^k 

(^ 


E(-^)l^; 


h.h I 


(81) 


As expected, this is the difference in the number of he¬ 
licity states. 

The physical number density and lepton density are 


1 


nh = 


"’41^ 


a(f)3 j (27r)3'^'^Al ’ 

dr’k 


1 f d^k 

( 82 ) 


The final lepton asymmetry is given by 


_ r '^L 

r]L = hm — 

t large 

TT 


1 




2C(3)r(f)3 a(f)3 y (2^)3 


E(-^)I4,J^ (83) 


which should be evaluated at a time after the comple¬ 
tion of reheating, so that a(t)T(t) approaches an asymp¬ 
totic constant value, but before electroweak sphalerons 
redistribute the charge between lepton and baryons. Fol¬ 
lowing this, further entropy production results in a final 
baryonic asymmetry about an order of magnitude smaller 
than r]L- 


VII. APPROXIMATIONS AND NUMERICAL 
ANALYSIS 

During the evolution of the Higgs VEV, c{r]) (defined 
in ([^)) is nonzero, which results in p{r]) 7 ^ 0 , signaling 
particle production. Additionally, since c{r], h = +) ^ 
c{r], h = —) generically, we expect a nonzero lepton asym¬ 
metry. At late times, 0 ( 77 ) —>■ 0 for both helicity values, 
resulting in /3 ^ exp(i|fe|? 7 ), which gives a nonzero asymp¬ 
totic value for |/3p, which is not generically identical for 
the two helicity values. Therefore, we expect a nonzero 
asymmetry to survive at late times, after the Higgs VEV 
(and hence Ml and ]jl) approaches zero. 

Calculating this asymmetry is complicated by the lack 
of analytic closed form solutions to the differential equa¬ 
tions (761, which must be solved numerically. In this 


section, we introduce a sequence of useful approxima¬ 
tions which simplify this problem significantly; we then 
present a numerical analysis of the resulting asymmetry. 
We focus particularly on the range of parameter space in 
which an asymmetry matching the observed cosmological 
baryonic abundance is generated. 


A. Higgs Oscillations 

We first note that it is desirable to have significant 
damping in the oscillations of the Higgs VEV. This is 
because the chemical potential poB vv' changes sign 
frequently during the oscillations, and so whether par¬ 
ticle or antiparticle production is favored also oscillates. 
Therefore, a significant damping in the amplitude of the 
oscillation avoids washout from this alternation. 
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As we explain below, this allows us to make two sim¬ 
plifications: First, that the asymmetry production occurs 
on a time scale during which a(t) is approximately con¬ 
stant, and second, particle production occurs primarily 
in those comoving momenta least affected by washout. 

We noted above that washout is significant unless the 
Higgs VEV is significantly damped. Consequently, the 
asymmetry production is dominated by the particle pro¬ 
duction during the initial relaxation of the Higgs VEV, 
which may be a fast process, compared to the evolution 
of the universe. We have mentioned in Sec. Ill Al that 
there are several reasons why the Higgs field may have 
a large vacuum expectation value after inflation. In one 
scenario, the Higgs VEV grows due to quantum fluctua¬ 
tions within the unmodihed Standard Model, or alterna¬ 
tively, the Higgs held may be trapped in a false vacuum 
during inhation. In the latter example, it is quite nat¬ 
ural that the evolution of the Higgs VEV, once it is re¬ 
leased from the false vacuum, would occur on time scales 
T <C \/H. This is more difficult to arrange in the former 
scenario, as the condition for the VEV to grow requires 
rueff < Hj, and the time scale of the Higgs VEV evolu¬ 
tion is ~ 1/meff. Rapid evolution of the Higgs VEV may 
still be arranged, as both iries and H are functions of 
time, although this may be somewhat unnatural. 

In the limit that the evolution of the Higgs VEV is 
rapid compared to the expansion of the universe, we may 
approximate a{t) ~ a{ts) constant, where we dehne t = 
ts to be the time at which the Higgs VEV begins rolling 
significantly. The comoving momentum during the epoch 
of particle production is fc = a{ts)p, and it is convenient 
to express P and a as functions of physical time t instead 
of conformal time 77 . They obey the differential equations 

^ = -C{t)p - iLoa, 

^ = C{t)a + iujp, (84) 


with 

C{t) = ^ = A [Mi (M(t)7reff + /ieff) 

-i\p\ + hpesf , 


(85) 


We remind our readers that the untilded HeS, Ml, and 
u! are the physical, and not comoving, quantities. For 
self consistency, we drop the terms proportional to H(t ), 
giving 


C{t) 


h 


MlPcS — (IpI + hpeff) Ml 


( 86 ) 


We emphasize that these expressions involve the physical, 
not comoving, momentum. However, when the evolution 
of the Higgs VEV is fast, these are related by the constant 
factor a(<s); we use this assumption to write 


VL = 


[ d^P 

2C(3)T(f)3a(t)3 i (27r)3 


7ra(fs)3 


y^i-h)\Pp,h{tE)\" 


(87) 


where tE is the effective end of particle production. 
We emphasize that our assumption is that a{t) is ap¬ 
proximately constant while the neutrino asymmetry is 
produced (for ts < t < tE), which allows us to use 
|fc| ^ o(fs')|p| in the integral of Eq. (83). Once the 
Higgs VEV relaxes to zero, no further asymmetry is pro¬ 
duced; however, the physical volume continues expanding 
~ a(f)3. This is responsible for the factor of a{t)^ in the 
denominator, which may be large; that is, this equation 
continues to hold even when a{t)/a{ts) 3> 1 , provided 
that a{tE) a(ts). 

Next we observe that if |/3p ^ 1 at all times, we can 
approximate d « 1 and the relevant differential equation 
is simply 


P = f C{t)dt. (88) 

Jo 

Sample plots of uj{t)'^C{t) (the factor cancels the l/w^ 
dependence in Eq. ( [^ ) and C{t) are shown in Fig.[^ We 
note that Ml/^-bS ~ ylM‘^v/Mj\jv ~ and 

y1/MM ~ 10“^° GeV“^ is fixed by the observed neutrino 
masses differences. Therefore, it is not surprising that for 
these parameters the typical scale of /ieff ~ yis about 12 
orders of magnitude larger than the typical scale of Ml- 
We proceed to describe the qualitative behavior of these 
plots. 

The C{t) panel (top) has a sequence of sharp peaks 
for h = — 1, of alternating sign, while these peaks are 
absent for h = +1. These are a consequence of the l/w^ 
factor in C{t). We note that during the initial pass of the 
Higgs field towards zero, y{t) is positive, and therefore, 
w ~ IpI + hy has a significant cancellation when \p\ ~ y 
for h = —1, while for h = +1 these factors always add. 
Additionally, for h = —1, |p| -|- hy changes sign at this 
peak. As the scale of y is larger than the scale of the left- 
handed neutrino mass Ml, C{t) changes sign “within” 
this peak. This is responsible for the very sharp positive 
and negative peaks. When evaluating this integral, these 
peaks cancel to a signihcant precision. 

To understand the behavior better, we consider 
uJ^C(t), which eliminates the sharp peaks. This plot is 
shown in the bottom panel of Fig. [T] We observe that in 
contrast to the C{t) plot, the h = +1 functions generally 
have smaller magnitudes than their h = —1 counterparts. 
This is because C{t) oc |p| -I- hy, and so the is a partial 
cancellation that affects w also affects the overall magni¬ 
tude of C{t) for h = —1. 

Additionally, for h = —1, the factors of (|p| -I- hy) 
changes sign (since y is negative), which affects the Ml- 
This factor is responsible for the various sign changes in 
the bottom plot of Fig. even though the Higgs VEV 
is decreasing during the entire time shown. We also note 
that the sharp spikes occur when y = 0, and so, momen¬ 
tarily, these plots are dominated by the Ml term, which 
(at this time) happens to be much larger in magnitude. 

As a simpler toy model, we consider a Higgs field which 










11 




FIG. 1. C{t) (top) and u!{t)‘^C{t) (bottom). The solid lines 
are for h = +1, while the dashed lines are for h = —1. These 
plots cover the initial pass of the Higgs VEV to zero. For 
concreteness, we have used the IC-2 scenario of Ref. EE] 
with the parameters A/ = 10^^ GeV and F/ = 10® GeV, 
along with M, the scale in Oq, equal to 10^^ GeV. This is 
about four orders of magnitude larger than the initial Higgs 
VEV. pmax is the maximum of |^(t)|. Units are appropriate 
powers of GeV. 


be independent, and furthermore, which may be chosen 
independently of any parameters describing inflation and 
reheating. We do note that in the scenario in which the 
initially large Higgs field VEV is produced via quantum 
fluctuations, vq will be determined by the scale of in¬ 
flation, although it will also be affected by any higher 
dimensional operators that influence the Higgs potential. 

We also made these further approximations: First, we 
assnmed ft Fb, such that 


Ml ~ l^e (1 - cos(2Ht)), 

^^e"^"*cos(2m). (91) 

Although we want significant damping, this is self- 
consistent, as we must have > F// in order for the 
Higgs VEV to nndergo oscillatory motion. 

Next, we observe that the cPp integral is dominated 
by momenta |p| ~ p-max, where /Umax is the maximum of 
\n{t)\ This is because these momentum values suffer the 
least washout during the subsequent oscillations of the 
Higgs VEV. Therefore, we approximate 


r] 


TT a{tsf mLx 
2C(3)T(t)3 a(f)3 (2^)3 


(92) 


obeys the equation of motion. 


d'^v 


3H 


dv 

dt 


+ m^v + Fh 


dv 

dt 


= 0 , 


(89) 


along with the boundary condition u(to) = I'D: v(to) = 
0. This is easier to analyze numerically, as opposed 
to considering the Higgs potential with running, and 
temperature-dependent, coupling constants. Further¬ 
more, the Standard Model Higgs field decays primar¬ 
ily through non-perturbative effects [H |3l]. For self- 
consistency, we again assume the Hubble friction term 
is negligible. This has the approximate solution 

v(t) = cos(f2(t — to)), (90) 


where ft = y/m^ — T^^. We have taken t = 0 as the time 
at which the Higgs starts oscillating, and so ts = 0. 

We note that a particularly interesting scenario is the 
case in which the potential V{v) = Xv'^/A] this is well- 
motivated by the fact that the Standard Model potential 
at large VEVs is dominated by this term. In this case, the 
term rn^v in the equation of motion would be replaced 
by Aeffn3, and in the solution the cosine function would 
instead be cn(uot/A^/"‘), where cn is a Jacobi sinusoidal 
function [32]. 

The Higgs potential is not known at large VEVs; there¬ 
fore, we will consider the parameters m, and Fh to 


B. Low Scale of Inflation 


For simplicity, we assume coherent oscillations begin 
instantly at the end of the inflationary epoch. We nor¬ 
malize the scale factor to one at the end of inflation, when 
the coherent oscillations of the inflaton start. In the fol¬ 
lowing computations, we approximate that the higgs os¬ 
cillations also start at this time. This is a good approx¬ 
imation in the case when the Higgs VEV is prevented 
from rolling by Hubble friction; when the Higgs is in¬ 
stead trapped in a false vacuum, this approximation will 
only be valid for sufficicently small barriers. 

Next we consider the factor of a{t)^T{t)^ in the denom¬ 
inator of T] given by Eq. (921. We emphasize that a(t)T{t) 


is the value approached at relatively late times, well into 
the radiation dominated epoch after Higgs relaxation has 
ended, but before the Standard Model degrees of free¬ 
dom have decoupled. This is completely determined by 
the two inflationary parameters, the inflationary scale A/ 
and the decay rate of the inflaton, Fj, where the Hubble 
parameter during inflation is 


Hi = 



(93) 


Since we do not fix a specific model of inflation, we 
take these to be independent parameters. We note that 
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FIG. 3. C(t) for uo = M = fl = 10^^ GeV and Th = 
10^^ GeV, which governs the differential equations for (5. 
Units are GeV“^ for time and GeV for C{t). 



C. Numerical Example 


FIG. 2. Contours of asymptotic late time aT; the scale factor 
is normalized to 1 at the time the inflaton starts oscillating. 
In the gray area, F/ > Hi, and there is no inflationary epoch. 


a{t)T{t) becomes constant once reheating has completed, 
and the asymptotic value is shown in Figure We can 
obtain the scaling of this factor as a function of the pa¬ 
rameters of the inflationary sector, under the assumption 
of instantaneous inflaton decay at H = F/ and of instan¬ 
taneous thermalization. This gives: 


As a numerical example, we used vq = M = = 

10^^ GeV and Th = 10^^ GeV, which give^max ~ 4.6 • 
10^^ GeV. A plot of C{t) is shown in Fig. the above- 
mentioned spikes are small for these parameters. 

P asymptotically approaches 0.0008 (for h = -hi) and 
0.006 (for h = —1). We have verified that ,5 < .1 at 
all times, so that our approximation in equation (881 is 
reasonable. The resulting asymmetry is 


lO^s GeV^ 
a(t)3T(t)3 ■ 


(96) 


A4/3 

a{t)T{t)cc i/s 1 / 6 , (94) 

Mpi i j 

or in terms of the reheat temperature, 

1 TrH .Qt-x 

(a(t)T(t))3 ^ Aj- ^ ^ 

where we recall that the scale factor a has been normal¬ 
ized to one at the end of inflation, which we take to be 
simultaneous with the beginning of Higgs relaxation. We 
see that a large baryon asymmetry is obtained for a low 
inflationary scale. For this reason we consider values in 
the bottom left corner of Fig. which are characterized 
by a relatively low A/. By consistency, Trh then needs 
to be below A/. 

We note that altering A/ and F/ may modify the evolu¬ 
tion of the Higgs VEV, particularly if finite temperature 
corrections to the Higgs potential are significant. In par¬ 
ticular, in scenarios in which the Higgs VEV begins in 
a false minimum which is destabilized by thermal fluc¬ 
tuations, there is a minimum reheat temperature which 
constrains A/ and F/. 


We must ensure that the Higgs energy density does not 
dominate the energy density of the universe, causing 
additional inflationary expansion, which requires A/ > 
GeV. For the minimum value, the inflationary Hub¬ 
ble parameter is then Hj = 2 • lO® GeV. If we take 
F/ = 10^ GeV, then at late times a{t)T{t) GeV. 

The resulting lepton asymmetry rjL « 10“^. 

We note that as m ~ H ^ Hj, a realistic implementa¬ 
tion of these parameters would likely have the Higgs VEV 
trapped in a false vacuum. The relaxation of the Higgs 
field would then commence after the start of coherent os¬ 
cillations, or ts > tend of inflation- The asymmetry is then 
enhanced by the factor a(ts)3/a(tend of inflation)^ > 1- 

This asymmetry will be diluted by a factor of 30 due 
to further entropy production, and it will be distributed 
between baryons and leptons by electroweak sphalerons. 
Therefore final baryonic asymmetry is about one or two 
orders of magnitude smaller. 

In summary, we have shown an explicit numerical ex¬ 
ample in which the asymmetry generated through neu¬ 
trino production is more than sufficient to explain the 
cosmological baryon abundance. We have also seen 
that production of a large asymmetry requires significant 
damping of the oscillations of the Higgs VEV, and also 
favors a low inflationary scale. 
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VIII. CONCLUSIONS 

In this work, we have explored particle production dur¬ 
ing an epoch of post-inflationary Higgs relaxation, with a 
particular emphasis on the production of a lepton asym¬ 
metry, which can be converted into the observed baryonic 
asymmetry. Unlike in previously considered models mi, 
the asymmetry considered here is produced via the evo¬ 
lution Higgs condensate directly, and does not involve 
interactions in the plasma produced by inflaton decay. 
Therefore, these models do not require a fast reheating, 
and in fact, we have shown a low reheating scale is desir¬ 
able. 

In particular, we have introduced a specific Oq opera¬ 
tor which involves only Standard Model fields (although 
extensions of the Standard Model may be necessary to 
produce this operator). This operator produces an effec¬ 
tive chemical potential for lepton number. 

We have solved the equations of motion exactly, in¬ 
cluding both this operator and a time-dependent Majo- 
rana mass. We then used a Bogoliubov transformation to 
relate the time-dependent creation and annihilation op¬ 
erators to the corresponding operators fixed at the time 
when the Higgs relaxation began. The resulting Bogoli¬ 
ubov coefhcients describe the rate of neutrino production 
during Higgs relaxation. From this, we calculated the re¬ 
sulting lepton asymmetry. 

After completing this formal analysis, we performed a 
numerical analysis, using a simplified model for the Higgs 
condensate evolution. This emphasized the importance 
of rapid condensate decay, which suppresses washout due 
to the oscillating sign of the effective chemical potential, 
and also the low reheat scale. We developed an approx¬ 
imation scheme that smooths out the sharp peaks that 
occur when |p| « We hnally illustrated a choice of 
parameters for which the resulting asymmetry is compa¬ 
rable to the observed value. 

Our scenario differs significantly from other scenarios 
of leptogenesis. In particular, the asymmetry can be 
generated for reheat temperatures well below the right- 
handed neutrino masses. This paves the way for a su¬ 
persymmetric generalization of the model in which the 
problem of gravitino overproduction does not arise. Fur¬ 
thermore, the final asymmetry is not tied to the param¬ 
eters of the neutrino mass matrix as in thermal leptoge¬ 
nesis, and a successful leptogenesis is possible even for 
the neutrino masses above 0.2 eV, in which case thermal 
leptogenesis is stymied by excessive washout [55] . 
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Appendix A: The Origin of the Og Operator 

In this appendix, we discuss methods of generating the 
Oq operator 


In the Standard Model in a flat static spacetime, the AB J 
anomaly allows the relation 

(A2) 


where A and B are that SUl( 2) and Uy(I) gauge helds, 
respectively, and Ug is the number of fermion genera¬ 
tions. The substitution of Eq. (A2) into Eq. (Al) is valid 


when the decay of electroweak sphalerons is fast, as com¬ 
pared to the Hubble parameter. Otherwise, the term AT 
involves the Chern-Simons number density, which is not 
changed by Higgs relaxation unless the phase of the Higgs 
VEV evolves. 

As to coupling these gauge fields to the Higgs field, we 
note that an effective term of precisely this form can be 
generated within the Standard Model, using quark loops 
and the CP-violating phase of the CKM matrix |5U [5S] . 
This term is small due to the small Yukawa couplings and 
small CP-violating phase. However, such a term can also 
be generated by heavier states with a different source of 
CP violation. The scale in the denominator may be the 
temperature, due to thermal loops, or the mass scale of 
new physics [55U57] . 

The sphaleron transition rate per unit volume at finite 
temperature, for constant Higgs VEVs, is 


Tsp = ka^T*exp{-Mw/gwT), 


(A3) 


where the exponential factor accounts for the suppres¬ 
sion due to being in the broken phase; it is equivalent 
to exp(—u/2T) where v is the Higgs VEV. Electroweak 
sphalerons are in equilibrium when this is greater than 
H^, where H is the Hubble parameter. The transition 
rate in the presence of a quickly evolving Higgs VEV 
has not been explored, although the rate during the elec¬ 
troweak phase transition from a = 0 to u = 247 GeV has 
been analyzed on the lattice, as a function of v{T) [58] . 

In section we found that the asymmetry is sup¬ 
pressed by a factor of {a{t)T{t))^, which favors a low 
inflationary scale. This generally corresponds to a slow 
reheating, while Higgs relaxation frequently occurs on 
a faster time scale. Therefore, during much of the re¬ 
laxation period, a > T and the sphalerons may not be 
in thermal equilibrium; the conditions for electroweak 
sphalerons to be in thermal equilibrium in the presence 
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of a time-dependent background have not been exten¬ 
sively explored. 

In each oscillation of the Higgs VEV, there is a brief pe¬ 
riod as the VEV passes zero during which v T, during 
which the above-mentioned suppression is absent. This 
also corresponds to the time of maximal particle produc¬ 
tion, which occurs when v = a{t)v{t) « 0. However, the 
time when the maximal asymmetry is produced is slightly 
offset from this time, as the effective chemical potential 
fl oc vv' is zero when u = 0. It seems unlikely that the 
time scale of sphaleron transitions will be less than the 
relevant time scale during which v T, even if the time 
of maximal asymmetry production is within this period. 
At the very least, it is difficult to arrange for this to hold. 

Therefore we note that, if there is another gauge group 
which couples chirally to leptons, it will also contribute 
to the divergence in equation ( |A2[ ). (The chiral coupling 
is necessary due to Furry’s Theorem.) Provided that 
interactions between the gauge field configurations and 
fermions are in thermal equilibrium, we find 

= (EW anomaly) + (A4) 

where F is the new gauge field and C is a constant deter¬ 
mined by the charges of the leptons and baryons under 
the new gauge group. Provided that these gauge bosons 
acquire masses which are not proportional to the Higgs 
VEV, it is possible for these to be in thermal equilib¬ 
rium at the relevant temperatures. (There may dynami¬ 
cal symmtry breaking in this sector, via a separate Higgs 
mechanism, or in the case of a t/(l) symmetry, via the 


Stiickleberg mechanism.) This equation can be rewritten 
as, 


ngCg^ 

327r2 




(AS) 


where jcs is the current associated with the elec- 
troweak Chern-Simons charge density. If the electroweak 
sphalerons are out of equilibrium, this is conserved, and 
therefor e has no effect on the analysis of sections |HI| 
through VH (similarly to how the baryonic current has 
no effect)) 

Therefore, if the electroweak sphaleron rate is insuffi¬ 
cient, we can couple the Higgs boson to a new gauge field 
combination, ea/ 3 /ii/F^'^F“^, to generate a term similar to 
(Al). As in the electroweak case, the coupling of 4)^ to 
can be accomplished through either thermal 
loops or heavy fermions. In the latter case, it is impor¬ 
tant that the fermions do not acquire masses through the 
Standard Model Higgs mechanism; otherwise, the Higgs 
VEV dependence cancels out. Such fermions may have 
soft masses similar to higgsinos and gauginos in super- 
symmetric models, or if a different Higgs sector is used 
to give masses to the F gauge boson, this field may also 
give masses to the relevant fermions. 

The divergence equation (A4) holds only in static, flat 


spacetime; the situation is more complicated in a curved 
and/or expanding spacetime. Generically, there may be 
contributions on the right hand side of the anomaly equa¬ 
tion, proportional to the gravitational anomaly |39j . 

If there are Nug right-handed neutrinos present, then 
(generalizing the results of [39] ') 


rigCg'^ 

327r2 




fj.jB+L-CS 


327r2 


(1-iV) 


a/375 ap-yS i i 

+ -^Sp.^Ss-c + + -(5“5^5a);. 


1 

6 



5a - 2 



(A6) 


where S describes the torsion of the spacetime, and R 
is the usual Ricci scalar. However, if there are the same 
number of right-handed and left-handed neutrinos then 


Appendix B: Conformal Higgs Field Equation of 
Motion 

Although we will use a toy model for our numerical 
analysis, it is beneficial to find the equation of motion for 
the comoving VEV v. From the Lagrangian in section [H] 


UgCg'^ 

32^2 


p/iiypa/3 


fijs+L-CS 


(A7) 


We consider only the scenario with = 1; that is, there 
are the same number of right-handed and left-handed 
neutrinos. 


dv 


= 0 , 


(Bl) 


where the derivatives signified with a prime are with re¬ 
spect to T] and V = a^V. We note that this is equivalent 
to the differential equation for the Higgs VEV v 


d?v 

dt'^ 


+ 3H 


^dv 

dt 


dV 

dv 


(B2) 
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It is necessary to express the potential in terms of comov¬ 
ing fields; as an example, we will do this with the 1-loop 
Standard Model Higgs potential, including finite temper¬ 
ature corrections. (However, in our numerical analysis. 


we will make use of a simpler effective potential for the 
evolution of the Higgs field.) The one-loop, zero temper¬ 
ature potential V, times a^, can be written 


K 


1—loop _ 




= +4^ 


(4- 


TT 


3a‘^mw(v)‘^ 


In 


mw{vf 




In 


ia^mzivY 


(mg^) _3_y 





Sa'^rritivY 






, (B3) 


r 


where S is the renormalization scale and the physical 
masses for the Higgs boson, Goldstone mode, W bosons, 
Z boson, and top masses are 


These have the same functional dependence on v as the 
regular masses have on v. Then the one-loop potential 
can be written: 


i-w 


2 2 
9 V 

4 ’ 

2 


2 {g^+g'Yv^ 

mz = - 


4 

^ -I-, mQ = rri^ + Xv^. 


mt = 


ytv 

V2’ 

(B4) 


It is convenient to define a comoving renormalization 
scale, S = aS, along with comoving masses 


~2 2 5^^^ 

mw = a —— 


2 ~2 
9 V 




ytv ytv 

""‘““72 “ V2 

rh^j = + 2Xv^) = + 2Xv^ 

rriQ = a^{m\ + Xv^) = mi + Xv'^. 


(B5) 


T/l-loop _ -2 , I ^ 

- 2™^"^ + 4^ + (47r)2 


mnivY ~ _ 3 

rht{vY 


^ (l„ (AfY - l) - (l. ^ 


(B6) 


We note that care must be used in evaluating the running couplings as functions of the comoving fields. During 
reheating finite temperature corrections may also be relevant; in terms of the comoving fields, these are 


Vt{v,T) = - 


2tt'^ 

f2 

2tt'^ 


6a mw{v) Jb 


2j f'mw{v)\ 2 r \2 T ('^z{v)\ , 2 . n2 T ('^t{v) 


V T 




+ 3a mz{v) Jb ^—j + 12a mt{v) Jp y 

+ Irk.ivfJs (Ifl) + 12Mvffy 


(B7) 


where 


CXD ^ 

JB{y) = V -^K2{ny), 

riZ 


n—1 


Jpiy) = ^- K2{ny). 


n—1 


and we have defined T = aT. Note the first three terms 


of Eq. (B8) and (B9) are fairly good approximation. 


(B8) 

(B9) 
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The relevant potential for the comoving VEV is then: and the projection operators are 


V{v, f) = + VT{i, f). (BIO) 




(C2) 


We note that as in EllS], it may be necessary to add fur¬ 
ther higher dimensional terms to the potential to produce 
a quasistable vacuum at large VEVs and/or to suppress 
isocurvature perturbations due to variations in baryon 
density. Additionally, dissipation effects may be relevant, 
and can also influence the production of a baryon asym¬ 
metry [40] , 


I note that 


1°Pl = 



7 “Ph =(";), (C3) 


and the complex conjugation operator is given by 


C = ^7^7° 






0 

-02 


• (C4) 


We also use the four-vector of Pauli matrices, 

a'^ = (l,a*) a^ = (l,-aO (C5) 

Appendix C: Two-Component Spinor Conventions 

Appendix D: Diagonalizing the Hamiltonian 

In the chiral basis, the Dirac 7 matrices are (Constant Mass and Chemical Potential) 


/ = 


(0 

A 

i 1 

^ 0 ctA 

1 

o)> 

7 = 1 

V-a* Oi 


In this appendix, we present the important steps lead- 
(Cl) ing from equation (55) to (56), for the interested reader. 
The two terms in (55) can be written as 


d^x v\ doVL = 


d?k 

(27r)3 


^ ^[u{h,k)*dou{h,k)] 


Ah)iih) (h)i (h) 




h,h 


-[vih, k)dou{h, 




(Dl) 


where we use the notation for the four-vector [E, —k). Similarly, 


d^x{dovl)i>L = I TbJa H ^)^ouih, k)*]af af'>- HKkD)*dou{h,k)*]af^'’af^'’x'i^'’x 

^ ’ h,h 

-[u(h, kD)dov{h, k)*dovih, k)]af^ 


-h) 

-k 


(D2) 


To evaluate the products of the spinors, we note that 


4 ^^^4 ^ ^-h-h = ^h,h, 

x(-'')(-fc) = C(fc,/i)x^"Hfc), (D3) 


where ^ is a phase that obeys: 

Ci-k,h) = -T]{k,h) 
a-k,-h)=v*Ck,h) 
ak,-h) = -fi*Ck,h). (D4) 

Additionally using the anticommutation relations, we 
may write the Hamiltonian as 
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H = - f [ [u{h, k)*dou{h, k) — u{h, k)dou{h, k)* — v{h, k)dov{h, k)* + v{h, k)*dov{h, 

2 j L k k 

— [u(/i, k)*dov{h, ko)* — v{h, kD)*dou{h, k)*]({k, — \v{h, k)dou{h, kn) — u{h, kD)dov{h, k)]C{k, 

K k,£) K]J K 

(D5) 


We note that the u’s and v’s depend on the momen¬ 
tum only through |fc|; therefore, u(h,ko) = u{h,k) and 
vih^ko) = v(h,k). The first combination is 

u{h, k)*dou{h, k) — u{h, k)dou{h, k)* + v{h, k)*dov{h, k) 
— v{h, k)dov{h, k)* = —2iCj [lap — |/3p] . (D6) 


while the other two are related by complex conjugation. 
One can show 


u{h, ku)*dov{h, k)* — v(h, k)*dou{h, ko)* = 2iCjha*P*. 

(D7) 


Together these give Eq. (56). 


Appendix E: Differential Equations for Bogoliubov 
Coefficients 


In this appendix, we present the important steps in de¬ 


riving the differential equations (64) from the equations 

(El) 


of motion. First, we introduce the notation 
9± = \/l ± / ■ 


The equations of motion require 

.du ,.-. _ , - 

i -—h n\k\u + UeffU = nAiLV, 
dr] 

dv^ ~ ~ 

i— -h h\k\v* + jlesv* = —hM^u*. 

dr] 

Since a, /?, and ui are time-dependent 


(E2) 


du a rv rjjfj /3 j urifi du da du d/3 a dn_ _ rv /3 dg+ rv 

dg -J2 \f2 da dg df3 dg y/2 dg -/2 dg 

dv .^ha -i n .-h/d i n dv da dv d(3 ha dg+ n h/d dg- , rv r.^^, 

dg -J2 \f2 da dg d/d dg y/2 dg ^/2 dg 

However, these functions also satisfy 

—iw-^g-e ~'‘-p iw-^g+e^ ^0 -I- h\k\u{h, k) -\- )lesu{h, k) = hMLv{h, k) 

V 2 v 2 

—ioj—^g+e~'^-^o — iui—^g-e^-^o -I- h\k\v{h, k)* + iiesv{h, k)* = —hMLu{h, k)*, 

V 2 V 2 

which allows us to simplify the above equations of motion to 

— — o_ e-* ■/'o" -t — ^ e* ■/'o" -f — ^ e-* ■/'o" -f ^ ^ e* J'o = q 
\/2dg^ ^/2dg^ ^/2 dg \/2 dg 

^—g^e-^dS^dn _ h^dd+^-if^Cbdr] _ u>df, ^ 


^/2 dg^^" \/2dg^ ' y/2 dg 

We may transform these into equations for the derivatives of a and fd 

1 da 
\/2 dg 


\f2 dg 


y It - T2 y w w) 


S didri 

? 

0 djdg 


(E3) 

(E4) 


(E5) 


(E6) 

(E7) 

(E8) 


3 If F is the antiderivative of dJ, then /J> uidg = F/v) - F{0). Differentiating this with respect to r, then gives dj(v). 
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We note that 


ff++5- — 1 + / + 1~/ — 2, 


and 


9- 


dg- 

drj 


9+ 


d9+ 

dr] 


1 + 

2 \ dr] dg) 


= 0 . 


The remaining combination is 


d9+ _ ^ d^ _ 1 df 

dg dg _ p dg' 


(Ell) 


(E9) 


Therefore, these equations simplify to 

da _ j3 1 df 2 i n Qdfj 
dg 2 ^1 - /2 dg 

dp ^ a I df fv Qdfj 

dg 2 - /2 dg 


(E12) 


(ElO) Since / = {h\k\ + /loff)/w, y/l - P = Ml/w, using the 
definition of uj above. Differentiating / gives us 


da 

dg 


P UJ 

'2i^ 


1 dfjQ^ h I I T /XefF ( / ~ \ dflefi dhdi^ 

-- [h\k\ + Heft)— -h Ml- 


dp a u) 
dg 2 Mj^ 


UJ dg uj^ 

1 d^efF /l I A; I ^efF 

UJ dg uj^ 


dg 


dg 
dML\ 


(h|fc|+/iefr)^+Mi — 
dg dg 


g2* JJ’ i^dri 




(E13) 


Finally, we can combine the dfies/dg terms, using — 
{h\p\ + pesY = ^L- T his gives Eq. (64), as desired. 


Applying equations (|E8|) in equations (E4) gives 


^ """ (E14) 

dg y/2 y/2 ^ ’ 

^ (E15) 

dg V2 v2 

which shows that the diagonalization of the Hamiltonian 
proceeds as in the time-independent case. 


I 

ticommutation relations, and the fact that uir^ku) = 
u{r,k) and v{r,ku) = v{r,k) because the three- 
momentum only appears as |fe| inside u and v allows us 
to write 


TefF = 


d^k 

( 271 -) 


3 


[d-p - i«n4>V’'' 




—uV a- '' a\ Cikjr) + vuaP'\P'\(k, r)* 
k ko k ko 


(El) 


Appendix F: Effective Lepton Number Operator 


In this appendix, we derive Eq. (0, starting from 
Eq. (70). Using the orthonormality of the spinors, an- 


Next we evaluate the products of the u’s and u’s, 


TefF — 


-h 


d^k 


E 


/ (|a|2 - |/3|2) ^ (^aP*e-^^ ^dn ff Qdf,^ 


—h 


2uj 

Ml 


— (+ faP 
2uj \ ) 


f{k, h) ' a'r 'a 


ko k 


_ L (^Q,*2g2i/J'wdr, _ 

2uj \ ) 




(F2) 
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where / is {h\k\ — fleg)/u! as above. Using the transformation equations, we recognize that this is 

d^k r „ ,(h)t .(h) 


LeS = 


(27 




ALeff, 


(F3) 


where 
ALeff = 


J (27ra)3 ^wLV J k k 

h 

-h- f^y^ih)^(h) _ / *2g2i/J'<idf; _ ^*2^-2i udfj\ ^(fc^/i)4fot^(fot 

2 \ / Kd 2 V / ^ 


(F4) 


Inverting the transformation equations gives 




/ a -h/?C 

U^*C a* 


\ kn / 


We note that |a|^ + 2|ap|/3p + |,0|^ = (|ap 4 
by the normalization condition, so this gives 


| 2\2 


)2 = 1 


(F5) 


Using this, we find 


ALeff = — 


E 


d^l Ml 


(2?™)^ uj 


L ^ k ko 


^ (27ra)^ 2w ^ 

g2i/J>i;d,; (_2|an/3|2 - |/3|4 _ |a|4) 

(-2|an/3|2 _ |a|4 - |/3|4) CAf'> 


_(_g-2i/J' ‘^d.rj j^{h) j^(h) 
k £) k 


(F7) 


Substituting this into Eq. (F3| gives Eq. 
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